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Abstract. Today’s high-intensity laser facilities produce short pulses can, in 
tight focus, reach peak intensities of 10^^ Wcm“^ and, in long focus, wakefield- 
accelerate electrons to GeV energies. The radiation-reaction-dominated regime, 
where the recoil from stochastic photon emission becomes significant, can be 
reached in the collision of such an electron beam with an intense short pulse. 
Measuring the total energy emitted in gamma rays or the presence of a prominent 
depletion zone in the electron beam’s post-collision energy spectrum would provide 
strong evidence of radiation reaction, provided enough electrons penetrate the 
region of highest laser intensity. Constraints on the accuracy of timing necessary 
to achieve this are given for a head-on collision. 
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1. Introduction 


The development of mnlti-PW laser systems, and the prospect of focussing laser pulses 
to intensities greater than 10^^ Wcm' mm, has led to substantial interest in how 
radiation reaction and strong-field quantum electrodynamics (QED) processes will 
alter the plasma physics studied over the next decade. Euture laser systems will be 
powerful enough such that the high fluxes of gamma rays (i.e. photons in the MeV 
range) and electron-positron pairs they produce in laser-laser [31 [H [5] and laser- 
solid m m interactions exert significant feedback on the evolution of the plasma, 
increasing the absorption of laser energy [14]. It will become possible to test the 
physics underpinning exotic astrophysical phenomena, e.g. pair cascades in pulsar 
magnetospheres mm, and illuminate fundamental questions about radiation reaction 
and the quantum vacuum in QED m- 

The importance of strong-held QED effects may be quantihed with the parameter 


rj = ^ ^ 

mcEsch 


( 1 ) 
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Figure 1. The experimental configuration under consideration: a GeV electron 
beam (blue) colliding with an laser pulse of intensity 10^^ Wcm“^ (red). The 
laser moves in the positive ^-direction; the electron beam centre is located in, and 
the electrons’ initial momenta are parallel to, the x-y plane. Observing radiation 
reaction will depend upon the degree to which the beams overlap, i.e. the collision 
timing and angle 0. The timing error is parameterised by (Ax, Ay), the location 
of the electron beam centre at t = 0, when the laser pulse is focussed at the 
origin. We do not consider the effects of a displacement between the beams in the 
z-direction. 


where is the electromagnetic field tensor and = jm{c, v) the electron four- 
momentum [3]. (This parameter is also called x m or Xe [ 12 ] •) ^ compares the 
magnitude of the electric field in the electron’s rest frame to that of the critical field 
of QED £^sch = 1.3 X 10^^ Vcm“^, which can produce electron-positron pairs directly 
from the vacuum m- 

T] also controls the stochasticity of photon emission. As the typical energy of a 
radiated photon for r] < 1 is huj OAArj^mc^ |3|, when 77 ^ 1, an electron can lose a 
substantial fraction of its energy in a single emission. Radiation reaction must be then 
treated as stochastic, rather than as the continuous loss of energy predicted classically. 
Pair production can become significant in laser-solid interactions laiT] for 77 ^ 1 ; in 
a colliding beams experiment, where the electron is not reaccelerated by the laser, 
77 ^ 1 is required for these processes to be significant [1210 • 

For highly relativistic electrons, 0 becomes 


77 - 


7|Ej 


■ V X B| 


-^Sch 


( 2 ) 


where 7 is the electron Lorentz factor, B the magnetic field and the electric 
field component perpendicular to v. In a laser-plasma interaction we would expect 
the electrons to have 7 ao, where oq = [/(A//im)^/1.37 x 10 ^^ is the 

laser’s strength parameter, / its intensity and A its wavelength. Then the typical 
77 ^ /A/(5.65 X 10^^ Wcm“^/im) and intensities > 10^^ Wcm“^ would be necessary 
to observe strong-field QED effects. 

However, if the electrons are pre-accelerated to high energies before encountering 
the high field region, the QED-dominated regime can be reached using much lower 
laser intensities. This was exploited in experiment E-144 at the SLAG facility [TSlfTB] . 
in which the collision of a 46.6 GeV electron beam and laser pulse of intensity 
10 ^^ Wcm“^ was observed to produce electron-positron pairs by photon-photon 
scattering. The advent of wakefield-accelerated electron beams with energies of 




Measuring quantum radiation reaction in laser-electron-heam collisions 


3 


1 GeV [H [H] or greater EOl HB and intense short pulses [26] now raises the 
possibility that this experiment could be repeated with an all-optical setup in today’s 
laser facilities. This has already been accomplished for target laser pulses with 
ao — 1 [22l[23]; recently Sarri et al m demonstrated production of multi-MeV gamma 
rays by non-linear Thomson scattering of a ^400 MeV electron beam. If high-Z solid 
rather than a laser pulse is used as a target, wakefield-accelerated electron beams can 
also be used to generate dense electron-positron plasmas by bremsstrahlung [25] . 

For collisions at uq — 1, the electron recoil is negligible. However, for a 
GeV electron beam and laser pulse of intensity 10 ^^ Wcm“^, 77 ^ 1 and quantum 
radiation reaction will dominate the dynamics of the electron beam. There is now 
a large body of work considering how strong-held QED effects can be observed 
in the spectra of the emitted gamma rays and the electron beam of such an 
experiment [271|28l|29l|30l|3ll[32l[33]. 

In previous work m, we considered the collision of a GeV electron beam with a 
laser pulse of intensity 10 ^^ to 3 x 10 ^^ Wcm“^; the beams collided at an angle of 180° 
and were timed to coincide perfectly at the laser focus. Quantum radiation reaction 
could be experimentally diagnosed by measuring the increased yield of high-energy 
gamma rays or the reduced energy of the electron beam after the collision. 

Both these signatures are sensitive to the overlap between the electron beam 
and laser pulse, which places constraints on the accuracy of collision timing required. 
In this article we will consider how the collision timing, and the angle between the 
beams, affects the total gamma ray energy, the electron beam’s final energy spectrum 
and therefore the viability of diagnosing quantum radiation reaction in an experiment 
using today’s high intensity laser facilities. 

Section describes the algorithm by which these interactions are simulated; 
section describes the chosen experimental parameters of GeV electron beam and 
10 ^^ Wcm“^, 30 fs laser pulse. Section presents simulated gamma ray spectra and 
shows that the yield is maximised in a head-on collision that occurs 80 /am from the 
laser focal spot. Section [^ presents simulated electron spectra and shows that, if the 
collision is accurate enough, radiation reaction leads to a substantial energy loss in 
that part of the electron beam the laser passes through. As measuring that loss can 
be a robust signature of radiation reaction, we constrain how accurate a collision must 
be for it to be sufficiently large. 

2. Method of simulation 

The collision of a GeV electron beam with a laser pulse of intensity 10^^ Wcm“^ 
is simulated with a single-particle, Monte-Carlo code that includes the electron’s 
oscillatory motion in the laser fields and the strong-held QED process of synchrotron 
photon emission. As pair production by the resulting gamma rays is likely to be 
negligible for these parameters m, we can neglect it here. 

Eollowing ISaiMlEZI, the differential rate of photon emission is 

^ V3a ri F{ri,x) 
dtdx 27 rTc 7 X 

X = h\F^yk^\/2mcEsch^ where the four-frequency = (cj/c, k), is the equivalent 
of 77 for the photon; the Compton time tq = hjmc? and the quantum synchrotron 
function, discussed in Appendix A, is dehned for 0 < x ^ 7/2- Integrating (§ over 
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that range gives the total rate of emission dN/dt oc r]h{r])/j, where the function 


,4) 

Jo X 

is a monotonically decreasing function of 77 that satisfies lim^^o h{r]) = bn/3. 

This rate is calculated in the Furry picture of QED [38], including an external, 
unquantised electromagnetic field. The basis states derived for this field may be used 
to calculate Feynman rules for quantised interactions IMlIlQlIlT]. 

Where the photon formation length is much smaller the characteristic length of 
the fields, that emission may be treated as pointlike and the fields quasistatically; thus 
the rate can be calculated in an equivalent system that has the same value of 77 . As the 
photon formation length Lph = A/ao [4], we may treat the laser fields as quasistatic 
provided uq ^ 1. In Erber [35], the equivalent system used to derive (§ is a static 
magnetic field in the limit 5 ^ 0 , 7 ^ oc. 

The emission algorithm follows the work of [SIE 11 MII 12 and is similar to that 
described in [32l[43]. Photon emission is a stochastic process and governed by Poisson 
statistics; the probability P of emission in a field with optical depth r is P = 1 — e~^. 
At the start of the simulation, and following each emission, the electron is assigned 
a ‘final’ optical depth rf = — log(l — P) for pseudorandom P G [0,1]. The electron’s 
optical depth r is integrated along its trajectory until rf is reached, according to 


dt Jq dMx ^ 


^/3a T] 
27 rTc 7 


h{r]). 


(5) 


At this point, the photon energy hw = 2xmc^ jr] is selected by solving 

r^Fiv.x') 


r h{r]) = 


f 


X' 


dx' 


(6) 


for X: where the pseudorandom number r is uniformly distributed in [0,1]. The 
electron then recoils antiparallel to its motion, with change of momentum huj/c. 
Energy conservation therefore requires that a small amount of energy is transferred 
from the laser fields during emission; as described in [5] , this leads to a fractional error 
in the electron energy of A 7/7 oc I /7 that may be safely neglected. While synchrotron 
radiation is directed forward into a cone of opening angle r\j 1/7 mi , we can assume 
that the photon and electron momenta are collinear, as the electron is always highly 
relativistic. 

Between emissions, the electron trajectory is determined classically by integrating 
the Lorentz force exerted by the laser fields (neglecting the beam’s space charge field, 
as described in section 


3.2) 


dr 

dt 

dt 


P 

7 m 

-e(E + V X B). 


(7) 

(8) 


This is implemented using a Boris push algorithm [45], so E and B are calculated at 
intervals staggered by half a timestep. 

The electron trajectories are discretised into timesteps of 0.005 fs, and at least 
10 ^ trajectories are followed for each set of parameters. 
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3. Parameters 

A schematic of the experimental configuration under consideration is given in figure 
We consider varying the angle of collision 0 as well as the mistiming in the x-y plane. 
This is parameterised by the offsets {Ax, Ay), which correspond to the displacement 
of the electron beam centre from the origin at t = 0 , at which time the laser is focussed 
there. 

3.1. Laser pulse 

The laser is a Gaussian focussed beam with waist wq = 2 jam and wavelength 
A = 0.8 jam, linearly polarised along the y axis and propagating towards -\-x. It 
has a Gaussian temporal intensity profile with maximum / = 10^^ Wcm“^ (ao = 65) 
and full length at half maximum / = c • 30 fs = 9 jam. The peak field Eq is related to 
the laser intensity by / = |ceoT’o- 

It can be shown from Q that rj, which controls photon emission, does not depend 
on the polarisation of the EM field, but only on the angle between the optical axis and 
the electron momentum 0. As such, the choice of polarisation here is arbitrary except 
that it defines the plane in which the electrons oscillate and so affects the angular 
distribution of radiation. This will be a fan of opening angle ^ao/y oriented along 
the polarisation axis [34]. 

3.2. Electron beam 

The energy distribution of the electron beam (up to a normalising constant) is 

for 100 MeV < £ < /i + cr/3; it has a peak at /r = 1000 MeV and a width a of 
250 MeV, which is related to the full width at half maximum (fwhm) / by / = 0.9cr. 
This is shown in figure]^ This distribution has been chosen as high energy wakefield- 
accelerated beams typically have large tails that extend to low energy. They also 
exhibit an angular divergence of order mrad, but here their initial momenta are 
directed along po = — (cos^, sin^, 0 ). 

The charge density of the beam forms a Gaussian envelope with (7\\ = 3.8 jam and 
cr_L = 4.2 jam (fwhm|| = 9 jam and fwhm^ = 10 jam) where the widths are given in the 
directions parallel and perpendicular to po respectively. 

The electric and magnetic components of the space charge force Esc — 
Qe/{'y‘^eoa^) of an electron bunch that has cylindrical symmetry nearly cancel if 
the electrons are ultrarelativistic 7 ^ 1. The force of radiation reaction E^r — 
ar^mcjTQ [3]. For Q 100 pC, a± 4 jam and 77 0.1, Esc ^ Frr means that 

the effect of the space charge field can be neglected and each electron’s trajectory in 
the laser pulse may be considered separately. 

Furthermore, as the typical separation between electrons {(j\\CF\e/Q)^^^ 5 nm 

(equivalent to 40 eV), we can also neglect collective radiation effects when considering 
gamma ray emission. 




Figure 2. The energy radiated by a GeV electron beam, per electron, to photons 
of all energies (left) and to photons with hjo > 500 MeV (right), when colliding 
with a Gaussian laser pulse of intensity 10^^ Wcm“^ at an angle 9 with parallel 
and perpendicular offsets Ax and Ay. 


4. Gamma ray production 

We consider the interaction of an intense laser pulse with a GeV electron beam, with 
parameters as described in section Figure shows the mean energy lost by such 
an electron beam to gamma rays when colliding with the specified laser pulse. 

We find that the total loss reaches its maximum of 120 MeV per electron for a 
head-on collision that occurs a distance Ax = 80 jam along the optical axis from the 
laser focus, but that the loss to photons with energy > 500 MeV is maximised at 
0.48 MeV per electron for a perfectly-timed collision that occurs at an angle of 50°. 

In determining the (Ax, Ay, 0) that maximise the energy emitted to gamma rays, 
we must consider three competing factors: a) the geometric factor in 77 ; b) each 
electron’s length of interaction with the laser pulse; and c) the overlap between the 
electron beam and laser pulse. 

a) arises from the definition of 77 (U): ifE and B are given by a linearly polarised 
plane wave and the electron is highly relativistic, travelling at angle 0 to the wave’s 
direction of propagation, then 

7 |E± + v X B| ^ 7(1 +cos 6 >)|E| 

■E’Sch -E'sch 

As 77 falls with increasing 0, colliding the electron with the laser pulse at an angle 
will have two effects: according to (§ the electron’s emission rate is reduced; and the 
photons it emits have lower energies, as the tail of the synchrotron spectrum grows 
non-linearly with increasing 77 . 

b) accounts for the importance of ‘straggling’ [5]: as photon emission is 
probabilistic, some electrons can reach the centre of the laser pulse having lost no 
energy; their 77 is then much higher and so too their probability of emitting a single 
photon with huu ^ ^m(?. In this case, the laser pulse is longer along the propagation 
axis than it is wide, so electrons traversing the pulse at an angle interact with the 
high intensity field over a shorter distance. While this reduces the total energy 
emitted, straggling becomes likelier with a shorter distance to reach the region of 
highest intensity; this should increase the yield of the highest energy photons. 

We can see the interplay of these two factors in the yield of the highest energy 
gamma rays, which is maximised for a collision at 50° rather than at 0°. This can be 
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Figure 3. The overlap between the laser pulse and electron beam (red, solid: 
left scale), normalised to the overlap at Ax = Ay = 0 = 0, and the total loss 
of energy to photons per electron (black, dashed: right scale), varying only Ax 
(left). Ay (centre) and 6 (right). 


explained by estimating analytically the average rj of an electron when it emits its first 
photon, which balances the reduced geometric factor against increased penetration into 
the laser pulse as 0 increases. The latter is likely to be more important than both 
the geometric factor and the overlap between the beams as it is the emission of the 
highest energy photons that is most sensitive to straggling [34]. The calculation is 
given in Appendix B; for the parameters under consideration here, 77 is maximised for 
a collision 0 = 65°, which is consistent with the peak shown in the right of figure]^ 
We conclude that it is most important to reduce the length of interaction and so the 
energy loss in the foot of the pulse to maximise the yield of the highest energy photons. 

However, maximising the total yield of gamma rays is achieved for a head-on 
collision at Ax = 80/im. We can understand this in terms of c), the overlap, which 
controls how many of the beam electrons interact with the laser pulse. Let us compare 
the total energy emitted in photons to an analytical estimate of the overlap 


n = 


hi 


d^rp{r,t)\E{r,t)\ 


( 11 ) 


where the beam charge density is as given in section 3.2 and |L^(r,t)| is as given in 
3.1, taking | sin0| = 2 / 7 r. That comparison is plotted in figure]^ varying each 


section 

of Ax, Ay and 0. 

We can see that the overlap models the dependence of the yield on Ax and Ay 
well. The yield increases with Ax up to Ax = 80 ym because the laser pulse diverges 
as it propagates away from its focal plane and therefore more electrons encounter the 
pulse. For Ax > 80 /im the yield falls because even though more electrons collide 
with the laser, they do so at lower peak intensity, which reduces their y and so their 
emissivity. 

Gamma ray production is more sensitive to the perpendicular displacement 
between the beams Ay than the parallel displacement Ax. This is because the peak 
intensity at the electron beam centre falls as exp(—2(A^/7Xo)^) for increasing Ay but 
as [1 -b (Ax/xr)^]“^ for increasing Ax; to reduce the peak intensity by a half requires 
Ax = 16 ym but only Ay = 1.4 ym. For Ay > 10 ym there is negligible photon 
emission because nearly all the electrons have missed the laser pulse. 

The overlap between the beams is nearly constant if only 0 is varied; however, the 
gamma ray yield is maximised at 80 MeV per electron at ^ = 35° and falls thereafter. 
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Figure 4. The energy spectrum of the electron beam before the collision (black, 
dashed) and after a collision at Ax = 0 (red) and 90 fim (blue). 


We can attribute this behaviour to the reduced length of interaction at large 6>, which 
reduces the energy loss before the electron reaches the pulse centre. However, this 
does not increase the gamma ray yield as much as increasing Ax. 

We conclude that allowing the laser pulse to diverge over a distance of 80 jam 
in a head-on collision is the optimal configuration to produce gamma rays. This can 
be achieved experimentally by focussing the high-intensity laser pulse with an optic 
that has an aperture in it. Provided this aperture is of sufficient size to permit the 
transmission of the wakefield-driving laser, the electron beam and resultant gamma 
rays, back-refiection and damage to the optical chain can be avoided. 

An analytical fit to the region in which the total gamma ray yield is at least 80% 
of its maximum value is 

Ax - 91 //m y / Ay V ^ ^ 

60/im / \3.6/amJ ~ 

for 0 = 0°. Even though the experiment would be designed for counterpropagation, 
the angle between the beams will vary from shot to shot due to the pointing variation 
of both the wakefield-driving and target laser pulses. Figure shows that to achieve 
any gamma ray production at Ax = 80 /im, that angle must be less than 15°. We can 
understand this requirement by considering the distance of closest approach between 
the centres of the electron beam and laser pulse Amin = Axsin6>/2; for Ax = 80/im, 
0 > 15° means Amin > 10 /im, the width of the electron beam, and so there is no 
overlap between the beams. 

5. Energy loss of the electron beam 

Detecting radiation reaction in an experiment could also be accomplished by 
comparing the energy spectrum of a wakefield-accelerated electron beam with and 
without the target short pulse, showing the reduced energy of the beam in the former 
case. Figure [^compares the energy distribution of the GeV electron beam described 
in section |3.2| to its initial energy spectrum for collisions at Ax = 0 and 90 /im. 

We can see in that in the former case the spectra are not sufficiently 
distinguishable, even though there is substantial loss of energy to gamma rays. That 
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Figure 5. The energy carried by the electron beam, per electron, per unit cross- 
sectional area a) prior to the collision and immediately after a collision at b) 
/\x = Ay = 0, c) /\x = 0, Ay = 5 ym and d) /\x = 50 ym. Ay = 0. 

loss is generated by those electrons that have collided with the intense part of the laser 
pulse; however, the peak at 1000 MeV remains as many electrons miss the laser pulse 
entirely. The broad energy spread is caused not only by the range of peak intensities 
encountered by the electron beam but by the stochastic nature of emission as well: 
two electrons travelling along the same trajectory will not necessarily lose the same 
energy. Neitz and Di Piazza [28] showed that this leads to energy broadening even for 
an electron beam without spatial extent and that radiation reaction manifests itself 
by ‘smearing out’ the initial energy distribution. 

If we consider the final state energy distribution for a collision at Ax = 90 /im, 
enough of the beam interacts with the laser such that the initial peak at 1000 MeV 
is entirely removed. If the electron beam is sufficiently well-characterised, it may 
be possible to detect radiation reaction by comparing the spectra from shots in the 
presence and absence of the target laser pulse. However, even high-quality wakefield- 
accelerated electron beams do not in general have energy spectra that are reproducible 
from shot to shot. An imperfectly generated electron beam would be a plausible origin 
of the blue spectrum in figure 

An alternative method of diagnosing radiation reaction exploits the fact that the 
laser pulse has a smaller diameter than the electron beam. Figure shows the energy 
carried by the electron beam over its cross-sectional area, i.e. its energy per unit 
volume integrated along its direction of propagation, for electrons that have collided 
with the laser pulse for various (Ax, A^, 0). We can see by comparing a) and b) that 
the laser pulse causes significant depletion of the energy spectrum in a region of radius 
2 /im around the optical axis. Resolving the electron’s beam areal energy spectrum 
could be accomplished by allowing to diverge over a long distance as it propagates 
away from the interaction region. A single shot would thereby allow the simultaneous 
measurement of the electron beam energy in the presence and absence of the target 
laser pulse and so the detection of radiation reaction. 

If we introduce a perpendicular displacement between the beams A// = 5 /im. 
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Figure 6. The reduction in the energy carried by the electron beam per unit 
cross-sectional area, at the point where the laser pulse passes through the beam. 


as in c), the depletion zone is still evident. However, d) shows that a longitudinal 
displacement, which allows the laser pulse to broaden before the collision takes place, 
makes the depletion zone less obvious. This is because the laser interacts with more 
electrons but at lower peak intensity, spreading the energy loss over a larger region of 
the electron beam. By considering the amount by which the energy is reduced at the 
centre of the depletion zone for collisions that take place at 6 > = 0 °, as in figurewe 
can constrain the accuracy of collision necessary for the depletion zone to be evident. 
An analytical fit to the region in which that loss > 3 MeV/rm“^ is 



(13) 


6. Conclusions 

The gamma ray yield and the energy loss of the electron beam are sensitive to both the 
maximal 77 reached by the electrons as well as the length over which this is sustained. 


However, for a realistic electron beam, the overlap between the beams will be a more 


significant factor in determining the strength of radiation reaction. 

Maximising these requires the counterpropagation of the electron beam and laser 
pulse; this exploits the slow angular divergence of the laser pulse, minimising the 
accuracy of longitudinal timing required, and the geometric factor (1 + cos^) in 77 . 

For a GeV electron beam colliding head-on with a 30 fs laser pulse of intensity 
10^^ Wcm“^, the mean energy lost to gamma rays will be maximised at 120 MeV 
per electron for a collision that occurs 80 jam along the optical axis from the laser 
focal plane. If the collision occurs closer to the focal spot, the gamma ray yield is 
reduced, but there will be a prominent depletion zone in the areal energy spectrum of 
the electron beam. 

The total energy emitted to gamma rays will be at least 80% of its maximum if 
the parallel and perpendicular offsets satisfy 



( 14 ) 
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and the depletion zone in the electron spectrum will be significant, i.e. the areal energy 
will be reduced by more than 3 MeV/im“^ from its unperturbed value, if 

2 

< 1. (15) 


Ay 


Ax 


4.6 /im \61 

The latter allows the simultaneous measurement of the electron beam energy in the 
presence and absence of the laser pulse and so the detection of radiation reaction. 

A head-on collision can be accomplished by focussing the high-intensity target 
laser pulse using an optic with an aperture of sufficient size to permit the transmission 
of the Wakefield-driving laser, the electron beam and resultant gamma rays. Obtaining 
a Ax and Ay that satisfy the accuracy requirements (14) and (15) will rely on gathering 
statistics over a large series of laser shots. While the region in Ax, Ay phase space is 
small, the energy converted to gamma rays will be sufficiently large, and the resulting 
depletion zone in the electron spectrum sufficiently prominent, to be robust signatures 
of radiation reaction for experimental parameters that can be achieved in current high 
intensity laser facilities. 
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Appendix A. The quantum synchrotron function 


The quantum synchrotron function is 

1 




— + 

71 1 - 2x/71 


POO 

K^/siS) - j 


dt 


where 



(A. 1 ) 


(A. 2 ) 


For small y, this is approximately 

F{v,x) = (f) ^ (A. 3 ) 

This means that ^( 77 , so the differential rate of photon emission (§, diverge 

as for low frequencies. However, the normalised energy emitted to photons with 

frequency between y and y + dy during interval dt, y^^dtdy, is well defined for all 
y. In particular, it is zero for y = 0. Similarly, the total rate of photon emission (§, 
which depends on /T(? 7 , y)/y dy, is always well defined. 


Appendix B. The typical 77 of an electron emitting its first photon 


Recalling that 


7|E_l + V X B| 

-^Sch 


7(1 -h cos6>)|E| 

^Sch 


(B.l) 
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Figure Bl. (left) The probability density that an electron colliding with 
Gaussian laser pulse of intensity 10^^ Wcm“^ at 0 = 30° (red) and 6 0° (b lue) 
emits its first photon at wave phase 4>i. The curves are calculated using ( |B.8| > and 
the vertical bars give the distribution found in simulations for a single electron 
with So = 1000 MeV colliding with a laser pulse as described in section [3T] The 
arrows indicate the value of ^i,mp given by ( |B.9| >. (right) Density plot of /(0i; 0): 
the red line is 0i,mp- 


for an ultrarelativistic electron propagating at angle ^ to a linearly polarised 
electromagnetic wave, the rate of change of optical depth against emission is 

dr 
dt 




(B.2) 


27rrc7 " 27rrc' £^sch 

The first approximation is that h{r]) ^ h(0) = Stt/S at all times along the electron 
trajectory. This is justified as in the foot of the pulse, the electric field is sufficiently 
small to keep rj < 0.1. Even then, h(O.l) = 0.93 h(0) so any overestimate of the rate 
will be negligible. 

In the presented simulations. So = 1000 MeV and therefore the electron 
can be treated as ultrarelativistic. Therefore its initial momentum is po = 
—fo/c (cos sin^, 0 ). As 70 = fo/^e ^ h can be assumed that any deflection 
of the electron by the laser pulse is minimal. Its trajectory before emission is linear, 
with position at time t given by {x^y^z) = —ct(cos 6 >, sin 6 >, 0). As (B. 2 ) depends on 
the electron energy only through ^( 7 ), which we have approximated as constant, the 
result will depend on So only in the sense that it must be sufficiently large that the 
electron be ultrarelativistic. 

In the simulations the full form of the Gaussian beam was used. However, for 
simplicity in these calculations, neglecting wavefront curvature and beam divergence 
gives 


E = Eo sin (/) exp 


21 n 2(^2 


(B.3) 


Wq 

(j) = kx — ujt = —(1 + cos 0)ujt (B.4) 

where uj is the laser’s angular frequency, k its wavevector, wo its waist size and / the 
full length at half maximum of its temporal intensity profile. 

Using (B.4) to recast (B. 2 ) and (B.3) in terms of 0, we find 
dr 


5a Eo 
d0 2y^UJTc ^Sch 


sm 




where 


d = 


■e 


2 In 2 


(1 + cos O^k'^WQ 


(B.5) 


(B. 6 ) 
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As the envelope in (B.5) varies slowly with we can approximate the oscillatory 
factor I sin 01 2/7r and thus integrate (B.5) from oo to 0i, the wave phase at which 

the electron emits its first photon. We find the following relation between the final 
optical depth and 0 i: 

_ ^ 5a Eq l-erf(C<?^i) 

2^/^U}Tc -E’sch C 

Photon emission is a stochastic process and governed by Poisson statistics, so 
these final optical depths are distributed as Tf ^ exp(—rf). Then the probability 
density that an electron colliding with a Gaussian laser pulse at angle 0 emits its first 
photon at phase 0 i is given by 

drf 


f{(t)i\0) = exp(-Tf) 


d0 


(B. 8 ) 


where we substitute (|B.7|) in the exponent and (B.5) in the modulus. Figure pT 


compares the analytical /( 0 i;^) with that obtained from simulation and finds good 
agreement. 

It is evident from figure |B 1 | that at larger angles of collision, the electron is 
more likely to penetrate further into the laser pulse before emitting its first photon; 
both the peak and the width of the distribution decrease as 0 grows, reducing the 
length of interaction. The most probable 0i = 0i,mp can be estimated by solving 
d/(0i; 6>)/d0i = 0, using (B.7), to find 

’ 1 / 5a 1 Eo y 

V VOtt wtc -Esch J 




1 


w 


(B.9) 


where W(x) is the Lambert W function, defined hy x = W{x)e^^^\ and is as given 

in (B. 6 ). 

We can substitute (B.9) into |E| ^ cxp(—and then r] = 7(1 + 
cos ^)|E|/£’sch to arrive at an estimate of the typical 77 of an electron emitting its 
first photon: 


T] - 


70(1 


i0)Eo 


-E’Sch 


exp(-C (/>i,mp) 


(B.IO) 


For a GeV electron colliding with the laser pulse described in section |3.H it is 
maximised for 0 = 65°. 
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